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We construct a scheme for the preparation, pairwise entanglement via exchange interaction, ma- 
nipulation, and measurement of individual group-II-like neutral atoms (Yb, Sr, etc.). Group-II-like 
atoms proffer important advantages over alkali metals, including long-lived optical-transition qubits 
that enable fast manipulation and measurement. Our scheme provides a promising approach for 
producing weighted graph states, entangled resources for quantum communication, and possible ap- 
plication to fundamental tests of Bell inequalities that close both detection and locality loopholes. 

PACS numbers: 03.67.Mn, 34.50.-s, 32.80.Wr, 03.65.Ud 



INTRODUCTION 



Entanglement is a vital resource for most quan- 
tum information processing (QIP) tasks, including long- 
distance quantum communication [2], teleportation- 
based quantum computation ^^^^ and one-way quantum 
computation (IWQC) [4]. An under-appreciated but cru- 
cial aspect of QIP is the need for speed of single qubit 
operations, to enable applications including synchroniza- 
tion of quantum communication networks, measurement 
and feed- forward in IWQC, and tests of local realism. 
For example, in IWQC, the processor speed primarily 
depends on the time needed for measurement and feed- 
forward, whereas the entanglement operation may be 
slow and accomplished simultaneously before commence- 
ment of the computation. In atomic systems, single- 
qubit fluorescence measurements are limited to microsec- 
onds due to auxiliary state lifetimes, and in alkali metals 
single-qubit rotation times are hampered by the gigahertz 
spectroscopic separations of hyperfine states. In this 
work, we overcome these obstacles by encoding in long- 
lived optical clock transitions (e.g., ^So ^ ^Po) of group- 
II-like neutral atoms, without sacrificing the advantages 
of other atomic schemes. Group Il-like atoms such as Yb 
and Sr have long been considered for atomic clocks and 
much recent experimental and theoretical effort has been 
dedicated to this group of atoms [i, @, 0, H, i, HO, El, El • 
The recent cooling of Yb into a Bose-Einstein conden- 
sate (BEC) and the ongoing study of interactions [12] 
make Yb an especially tantalizing candidate for atomic 
qubits. Our approach for entanglement and measure- 
ment of group-II atoms offers promising techniques for 
the high-speed synchronization needed for quantum com- 
munication and computing, and also for the near-term vi- 
olation of a Bell inequality in a single laboratory, without 



any assumptions about signaling, sampling, or enhance- 
ment [13, ^,[151, [3. 

Significant experimental progress has been achieved 
towards entangling atoms in optical lattices [17], which 
could lead to the creation of an initial state for IWQC. 
Here we take a complementary approach, considering the 
entanglement of individual pairs of atoms on demand, 
comparable to other addressable neutral atom architec- 
tures [18|. Rather than creating a generic cluster state, 
we propose the creation of computation-tailored weighted 
graph states as a resource for IWQC and other QIP 
tasks. Our technique combines efforts to prepare individ- 
ual atomic qubits from a BEC [19], coherently manipu- 
late and transport atoms [2O, HJusing optical tweezers 
at a "magic wavelength," entangle atoms via an inher- 
ently robust exchange interaction [22, 23], rotate single 
qubits via a three-photon optical dipole transition |1Q| |, 
and perform fast (~ns) measurements via resonantly en- 
hanced multi-photon ionization (REMPI). A "loop-hole 
free" Bell inequality test imposes stringent requirements 
on detector separatio n il5J a nd efficiency (see, e.g. the 
experimental work in [1^, [2J), and presents an enticing 
test-bed for fast measurements with applications to QIP. 
We study the limits of fast measurement for encoding in 
the optical clock states of Yb and Sr, which can be re- 
solved spectroscopically and measured on a ~10ns time 
scale, thereby admitting space-like separation over a few 
meters (as opposed to large spatial separations consid- 
ered in [24]). We show that such Bell tests in a single 
laboratory should be feasible via a detailed theoretical 
analysis accompanied by comprehensive numerical simu- 
lations. 



II. QUBIT PREPARATION AND TRANSPORT 



Clock transitions in ions have been considered for effec- 
tively encoding qubits for ion trap-quantum computing 
*Electronic address: restock@physics.utoroii to.ca| due to extremely low decoherence rates [26|. Simi- 
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larly, in the case of neutral atoms, optical clock transi- 
tions in alkaline- earth and group-II-like atoms are ap- 
pealing candidates for encoding qubits. Single atoms 
have been experimentally isolated [l9] and transported 
in optical dipole traps [2Q,, 21] • By trapping at a "magic 
wavelength" 0, [1], the light shift potential is made ef- 
fectively state-independent, ensuring phase stability of 
the qubits for several seconds. For example, for the clock 
states of Sr, the light shift dependencies on the trap laser 
frequency u differ by dA/du = 2.3 x 10"^° [t^. Therefore 
light shift fluctuations can be kept to less than 0.1 Hz by 
using a trap laser with line width of 100 MHz. Further- 
more, the magic wavelength at 813.5 nm (easily accessible 
using commercial lasers) is far detuned from the excited 
states so that photon scattering rates are on the order 
of 10 s for trap light intensities of 10 kW/cm^ [7]. This 
ensures a coherence time of 10 s or more for trapping and 
transporting atoms. 



III. ENTANGLING OPERATION 

We devise a universal entangling operation for bosons, 
analogous to the recently proposed fermionic spin- 
exchange gate [i^. This gate is based on the exchange 
interaction recently demonstrated for bosonic Rb atoms 
in a double- well optical lattice Because of inher- 

ent symmetrization requirements, gates based on this ex- 
change interaction offer a natural resistance to errors and 
greater flexibility for encoding atoms, thereby enabling 
an entangling operation even for atoms with interaction 
strengths that are state-independent (e.g., Rb [23]) or 
partially unknown, as is the case for most group-H-like 
atoms (e.g., Yb [l2l|). 

The entangling operation is achieved by temporarily 
bringing together a pair of atomic qubits via mobile 
optical tweezers. Unlike state-dependent optical traps 
wherein atoms are trivially separated into opposite wells 
after interaction, we have state-independent traps in 
which the dynamics of the system generally determine 
the likelihood of a successful separation. However, under 
adiabatic conditions the atoms definitely end up in oppo- 
site wells. We assume a strong confinement to one dimen- 
sion (ID) by higher order Hermite Gaussian beams ac- 
cording to [i^ . All the essential physics is captured in the 
ID model we employ here, although performance could 
conceivably be enhanced by exploiting multi-dimensional 
effects such as trap-induced resonances p7| . 

The Hamiltonian for two trapped atoms a and b with 
internal structure , G {|0), |1)}) is given by 

^= ^[Ha^Hb^2aijhu;^S{xa-Xb)]^\ij){ij\ (1) 

for Ha,b = Plb/2m + V{xa,b - d/2) + V{xa,b + d/2), 
with Xa,b and Pa^b the position and momentum of atom a 
or b. The tweezer potential V{x) = —Voexp{—x'^/2a'^) 
describes a Gaussian trap of depth Vq and variance a^. 



The two wells are separated by a distance d, u;± is 
the harmonic oscillation frequency of the transverse con- 
finement [2^, and aij is the state-dependent scattering 
length for the two-qubit states = \i)a^\j)b' 
merically solve the Hamiltonian dynamics of individual 
qubit states using a split-operator method. Two-atom 
energy spectra are plotted as a function of well separa- 
tion (Fig. [1]) for different interaction strengths. 

Due to symmetrization requirements, not all com- 
binations of vibrational and qubit states are allowed. 
For example, a pair of composite bosons cannot share 
the ground state if the qubits are in the antisymmet- 
ric state 1^-), defining = (|01) ± \10))V2. As 
in the fermionic case [22|, it is possible to exploit these 
symmetrization requirements in order to produce a two- 
qubit entangling operation for bosonic atoms (see [2£] 
for details). Consider a pair of identical bosons, one lo- 
calized in the left trap (|V^^)) and carrying a qubit in 
the state \(p^) = a|0) + /3|1), the other in the right 
trap (IV^^)) and carrying a qubit in the state \(p^) = 
/i|0) + z^|l). The initial symmetrized wavefunction (as 
a tensor product of vibrational and qubit states) is then 




012345 012345 



cy(units of a) c/(units of a) 

FIG. 1: (Colour online) Adiabatic energy levels as a function 
of well separation. Energies are measured in units of huo, 
where ujq is the harmonic oscillation frequency of one atom 
in a single well, (a) Lowest six energy levels for aij = O.lcr. 
Energy levels correspond to symmetric or antisymmetric ex- 
ternal eigenstates. The antisymmetric curves (red) are the 
lower of the two curves at ^ ^ — 7 hooo and the lowest of the 
three curves at E ^ —5.8 huuo for d = 0. (b) Lowest two levels 
of (a) for different scattering lengths. The lowest three energy 
curves (from bottom to top) correspond to aij = 0, a^ = O.lcr, 
and aij = cr, and asymptote to the antisymmetric (topmost) 
curve for infinite aij. The antisymmetric eigenstates are not 
affected by the interaction and hence the topmost (red) curve 
does not shift for different aij. 

As the wells are brought together and separated adi- 
abatically, the energies evolve as shown in Fig. (TJ and 
each two-qubit state |00), |11), and |^^) acquires a phase 
^00, </>ii5 and depending on its respective energy 
curve. Adiabaticity can be satisfied even for negative 
scattering lengths, since transitions between vibrational 
states of different symmetry or parity are suppressed. 
For constant tweezer speed the adiabaticity criterion 
is <C (jJiLo'^^/Vo. Here, hcOab is the energy difference 



3 



between any coupled states. Time-dependent numerical 
simulations confirm the validity of the adiabatic approx- 
imation over a wide range of values of Vq and aij [29^ . 
The final state after an adiabatic change of separation is 

l^f) = l^-)®(^e-^-|*-» (2) 
+ |V'+)®(aAie-''^™|00)+/3i/e-*'^"|ll) + ^^^e-^<^+|*+)), 

using |^/;±) = {\tp^i^^) ± \i^^i^^))/V2. 

Evidently this process corresponds to a tensor prod- 
uct of the identity acting on the vibrational state and a 
unitary U acting on the qubit state. Thus, the internal 
qubit evolution simplifies to 

U = e-'-^^lOOXOOl + ^ '^^^^ '^" (|01X01| + |10)(10|) 
+ ^^^^^Y^(|01){10| + |10)(01|) + e-^-|llXll|. (3) 

As in [s^, a controlled-phase gate can be obtained 
even if ^+ 7^ ^_ by sandwiching a single-qubit phase 
gate between a pair of U operations. That is, G = 
U [S{7t) 5'(0)] U for S{0) = exp(i(9|l)(l|). Thus defined, 
G is locally equivalent to exp(— i7|ll)(ll|) if 

^00 + -(/>+-(/>- = (2n±^)7, V n G Z. (4) 

As shown in Eq. ([2]), the phases critical to this en- 
tangling operation are acquired in a non-separable basis. 
This leads to the inherent robustness observed in initial 
experiments [23]. In standard schemes, the important 
non-separable phase is usually acquired due to the inter- 
nal state dependence of the interaction strengths aij. In 
the case of this exchange symmetry-based gate, however, 
there always is an energy gap between symmetric and 
antisymmetric curves. The singlet state |^~) therefore 
acquires a phase different from the triplet states even if 
the interaction strengths are state-independent (except 
as aij ±00). This substantial phase difference en- 
ables the exchange gate to operate faster than standard 
collisional gates that rely on the difference in aij. Fur- 
thermore, this gate works over a large range of scatter- 
ing lengths [see Fig[T](b)], which is especially important 
when designing experiments for atomic species with any 
currently unknown scattering lengths (e.g., Yb or Sr). 
Current studies of Yb interactions [l2[ already promise 
a wide applicability of this entanglement gate for differ- 
ent isotopes. (For ^^^Yb, aoo ^ 13 nm and for ^'^^Yb, 
aoo ^ 5.6 nm. aoi and an, are not yet known.) 

IV. SINGLE QUBIT ROTATION AND 
MEASUREMENT 

Recent attempts to cool and trap neutral Yb and Sr 
have been very successful, and we therefore consider them 
primarily. Optical clock states in Yb and Sr have ex- 
tremely low decoherence rates, due to the fact that elec- 
tric dipole one- and two-photon transitions between ^Sq 



and ^Pq states are dipole and parity- forbidden, respec- 
tively [see FigsISfa) and [3l^a)for energy levels and tran- 
sition wavelengths]. While affording long lifetimes, the 
selection rules also present a significant challenge to fast 
coherent manipulation and measurement of qubit s. To 
overcome this challenge, we employ a coherent, three- 
photon transition to perform single qubit operations, uti- 
lizing the excited ^Si and ^Pi states [10]. The three 
transitions ^Sq ^Pi, ^Pi ^Si, and ^Pq are 

electric-dipole allowed (see [sl, Q for transition matrix ele- 
ments). Because three beams can always be arranged in a 
plane such that the transferred recoil cancels, this three- 
photon transition has the benefit of being recoil- free [10|]. 
For Sr, the need for three lasers may be reduced to two, 
as explained below. 

We model this three-photon transition by a master 
equation using the Liouvillian matrix given in [10]. Its fi- 
delity is limited by the short-lived intermediate ^Si state, 
which decays primarily to the ^Pi state. The fast coher- 
ent rotation of qubits is followed by the fast readout of the 
^Pq state via REMPI on a nanosecond or even picosecond 
time scale. Re-using the excited state, photoioniza- 
tion can then be accomplished in a two-step process. An 
on-resonant ^Pq to ^Si transition is followed by a final 
ionization step at A < 563 nm for Yb and A < 592 nm 
for Sr. The main errors in this read-out scheme are due 
to population in the ^Pi to ^Si states. During readout, 
any population in ^Pq and ^Si will be counted as logical 
|1) (ionized). Population in and ^Pi will be counted 
as logical |0) (not ionized). 

The case of Sr is particularly interesting: the tran- 
sitions ^So Pi and ^Pi Si are close in energy 
difference (689 and 688 nm, respectively) so that a reso- 
nant two-photon transition ^5*0 -^^ Si utilizing a single 
laser is possible. This reduces the laser requirement from 
three to two. Figure [2](b) shows fidelities for qubit rota- 
tion for wavelengths in the range 688 to 689.5 nm. The 
time for a 7r-rotation is minimized by tuning to 688.7 
nm. Figure ^c) shows the fidelity and time scales for a 
TT-rotation as a function of laser powers. For fairly realis- 
tic mode- locked laser powers, 10^ W/cm^ (roughly 1 kW 
pulse peak power focused onto 100/im^), rotations within 
a few nanoseconds are possible with better than 90% fi- 
delity. Higher fidelities of 99.99% can be reached for the 
same detuning by using lower laser powers of 10^ W/cm^. 

V. TESTS OF LOCAL REALISM 

We show the efficacy of our fast measurement scheme 
by applying it to a test of local realism. This 
is expressed in the usual Clauser-Horne-Shimony-Holt 
(CHSH) form [3, 

(B) = (QS) + (RS) + (RT) - (QT) < 2, (5) 

for local realistic theories, whereas Tsirelson's quantum 
upper bound is 2V^. For a |^^) entangled state, the 
quantum bound is saturated for Q = R = X, 
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FIG. 2: (Colour online) (a) Energy levels of Sr and three- 
photon transition for manipulation of the qubit encoded in 
^5*0 and ^Po- (b) Minimum fidelity Fmin of single qubit oper- 
ation in Sr (solid red line) and time scale for tt pulse (dashed 
blue line) as a function of Ai (= A3) using a peak laser pulse ir- 
radiance of 10^ W/cm^. A2 is determined by the on-resonance 
condition for the three-photon transition, (c) Minimum fi- 
delity Fmin of single qubit operation (solid red line) and time 
scale for tt pulse (dashed blue line) as a function of laser irradi- 
ance /peak- Detuning is fixed to Ai = 688.7 nm. (d) Resulting 
expectation value of the Bell operator (B) and threshold for 
a local hidden variable model (solid black line). 



S = {X - Z)/V2, and T = (X + Z)/V2, with X, F, Z 
the Pauli operators. These measurements are obtained 
via basis rotations 7l{0) = exp(+z^Y'/2) applied to the 
state, followed by measurements in the z-basis. This 
corresponds to measurements of the form Q = UqZUq 
with Uq = 1, Ur = 7^(7^/2), Us = '^(37r/4), and 
[/t = 7^(7^/4). 



gled Bell state via a controlled phase gate as discussed 
above, then separating the atoms by a few meters. In 
a far-off-resonance, magic- wavelength trap, qubit coher- 
ence times are on the order of 10 s or longer. For accel- 
erations of 200 mm/s^ or faster [21], separations of a few 
meters should be feasible. At this distance, synchronous 
measurements on a nanosecond time scale are required 
to ensure space-like separation. Within this time win- 
dow, the measurement basis is chosen randomly, qubits 
are rotated to reflect the choice of measurement basis, 
and qubit states are measured in the computational ba- 
sis using REMPI. Fast random basis selection can be 
accomplished by using a light emitting diode (LED) as 
in [15]. The time necessary for this random basis selec- 
tion can be minimized [e.g., by using shorter signal paths 
and custom-built electro-optic modulators (EOMs)] to 
ensure basis selection times of less than 10 ns. Rotation 
of the measurement basis is achieved via a coherent cou- 
pling of the qubit states via three-photon Raman tran- 
sitions. The presence of the ion (i.e., the freed electron) 
will be detected via a single channel electron multiplier 
with above 99% efficiency [31] . 

As in a typical single channel experiment [16], the 
measurement outcome can be only "ion"= |1) or "no 
ion"= |0). No data are discarded, and no assumptions are 
made about "fair sampling" [13] or "enhancement" [l3|. 
Loss of an atom will result in a "no ion"= |0) count, 
which reduces the degree of Bell inequality violation but 
does not open any loopholes. High transport and detec- 
tor efficiencies are necessary to ensure that a violation 
occurs. A calculation of the CHSH-type Bell inequality 
violation [l3|, including errors in rotation and ionization 
readout, is shown in Fig. [SJ^d) for Sr and Fig. [31(b) for 
Yb. To achieve an average value of the Bell-operator 
larger than 2, as required for a violation, measurements 
on a time scale of a few nanoseconds (including signal 
processing times) should be possible with either atomic 
species. 



(a) (b) 




FIG. 3: (Colour online) (a) Energy levels of Yb and three- 
photon transition for manipulation of the qubit encoded in 
^5*0 and ^Pq. (b) Expectation value of the Bell operator for 
imperfect single-qubit rotations in Yb as a function of time 
scale of the measurement for /peak = 10^ W/cm^. 

Inequality ([5|) is tested by first preparing an entan- 



VL CONCLUSIONS 

We propose schemes for fast recoil-free manipulation 
and measurement of qubits in Sr or Yb, and discuss an 
entangling operation for identical bosons in optical tweez- 
ers based on the exchange interaction first discussed for 
fermions in [22]. We furthermore show that it is possible 
to simultaneously close both space-like separation and 
detection loopholes for group-II-like atomic qubits sepa- 
rated on only a laboratory scale. This lays the ground- 
work for future exploration of measurement-based com- 
putation. Finally, our work identifies major challenges 
and provides concrete guidelines for experiments utiliz- 
ing bosonic Yb or Sr for quantum information processing 
applications. 
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